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Abstract 

The aim of the Swedish-Japanese EISCAT campaign in February 1999 was to mea- 
sure the ionospheric parameters inside and outside the auroral arcs. The ion line 
radar experiment was optimised to probe the E-region and lower F-region with as 
high a speed as possible. Two extra channels were used for the plasma line measure- 
ments covering the same altitudes, giving a total of 3 upshifted and 3 downshifted 
frequency bands of 25 kHz each. For most of the time the shifted channels were 
tuned to 3 (both), 4 (up), 5.5 (down) and 6.5 (both) MHz. 

Weak plasma line signals are seen whenever the radar is probing the diffuse 
aurora, corresponding to the relatively low plasma frequencies. At times when 
auroral arcs pass the radar beam, significant increases in return power are observed. 
Many cases with simultaneously up and down shifted plasma lines are recorded. In 
spite of the rather active environment, the highly optimised measurements enable 
investigation of the properties of the plasma lines. 

A modified theoretical incoherent scatter spectrum is used to explain the mea- 
surements. The general trend is an upgoing field-aligned suprathermal current in 
the diffuse aurora, There are also cases with strong suprathermal currents indicated 
by large differences in signal strength between up- and downshifted plasma lines. 
A full fit of the combined ion and plasma line spektra resulted in suprathermal 
electron distributions consistent with models. 



1 Introduction 

The incoherent scatter spectrum consists mainly of two lines, the widely used and 
relatively strong ion line and the very weak easily forgotten broadband electron 
line. There is also another line present, the plasma line, due to scattering from 
high frequency electron waves, namely Langmuir waves. From the downgoing 
and upgoing Langmuir waves, two plasma lines can be detected by the radar. 
The frequency shift from the transmitted signal is the frequency of the scattered 
Langmuir wave plus the Doppler shift caused by electron drift. Plasma lines can 
be used to measure the electron drift and hence the line-of-sight electric current. 
The problem in ion line analysis with the uncertainty of the radar constant can be 
solved by the plasma line frequency determination and when that is done the speed 
of measurement can be significantly increased by including the plasma line in the 
ion line analysis. However, since the frequency of the plasma lines is not known 
beforehand, and the frequency is varying with height, it is difficult to measure 
them with enough resolution. 

There have been a number of reports on plasma line measurements and their 
interpretation. Most of them have discussed the frequency shift from the transmit- 



ted pulse and the scattering has mainly been from the F-region peak, e.g. Showen 



{1979)), [Kofman ct al. (1993| ) and [Nilsson ct al. (1996a| ). The latter two showed 



also that the simple formula for the Langmuir wave frequency, 

/ 2 = / 2 (l + 3fc 2 A 2 ) ) + / 2 sin 2 a (1) 

where f p is the plasma frequency, k the wave number f c the electron gyro frequency, 
\d the electron Debye length and a the angle between the scattering wave and 
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the magnetic field, is valid to within a few kHz and thus enough to set the radar 
system constant. To be able to deduce any electron drifts, or current, out of the 
positions of the lines these authors also show that E q. (|l|) is not sufficient, and it 
is necessary to carry out more accurate calculations. Hagfors and Lchtinen (1981 ) 
had also to go to further expansions in deriving the ambient electron temperature 
from the plasma lines. The fact that so many reports deal with the F region peak 
is due to the altitude profile shape of the plasma line frequency, which according 
to Eq. (|l|) will also show a peak around that height. The measurements can thus 
be made relatively easily using rather coarse height resolution but good frequency 
resolution and only detect the peak frequency. Measurements using the same 
strategy, but at other heights, have been made with a chirped radar by matching 
the plasma line frequency height gradient and the transmitter frequency gradient 
(Birkmaycr and Hagfors, 1986; [sham and Hagfors, 1993). This technique allows 
determination of the Langmuir frequency with very high frequency resolution, but 
do not use the radar optimally, since the chirped pulse cannot be used for anything 
else than plasma line measurements. 

The enhancement of the plasma lines, which occurs in the presence of suprather- 
mal electron fluxes ( Perkins and Salpeter, 1965| ), either p hotoelectrons or secon - 
daries from auroral electrons, has been investigated by Nilsson ct al. (1996b| ), 
where they also calculate predictions of plasma line strength for different incoher- 
ent scatter radars and altitudes. They also show that the power of the plasma line 
is rather structured with respect to ambient electron density, depending on fine 
structure in the suprathermal distributions due to excitations of different atmo- 
spheric constituents. 

Incoherent scatter plasma lines in aurora are more difficult to measure since 
the variations in the plasma parameters are strong with large time and spatial 
gradients. Reports of auroral plasma lines in the aurora are also more rare, and 



most of them are based on too coarse time resolution ( 


Wickwar, 1978; 


Kofman 


and Wickwar, 1980; 


Oran et al., 1981; 


Valladares et al., 1988 


), with resolutions 



ranging from 30 seconds up to 20 minutes. The enhancements over the thermal 
level were high, but consistent with what could be expected of model calculations 
of suprathermal electron flux. They also tried to calculate currents and electron 
temperature from the frequency shifts of the plasma lines but with very large 
error bars. Kirkwood et al. (1995) used the EISCAT radar and the filter bank 
technique and recorded much higher intensities of the plasma lines, since they got 
down to resolutions of 10 seconds, and showed also that the plasma-turbulence 
model proposed by Mishin and Schlegel (1994) was not consistent with the data, 
but could be explained by reasonable fluxes of suprathermal electrons. 

In this paper we present data obtained with the high resolution alternating code 
techniq ue (Lchtinen and Haggstrom, 1987 ), as was also done for F-region plasma 
lines by Guio ct al. (1996[) , with even higher intensities due to the time resolutions 
of 5 seconds. An interesting, but at the time of the experiment not realis able at 
EISCAT, techniq ue would have been the type of coded long pulses used by Sulzer 
and Fcjcr (1994 ) for HF-induced plasma lines. From relative strengths between 



up- and downshifted lines we detect a general trend of upgoing field-aligned cur- 
rents in the diffuse aurora carried by the suprathermal electrons. We propose a 
generalisation of the theoretical incoherent scatter spectrum, to include multiple 
shifted electron distributions, and in one example we do a full 7-parameter fit of 
the incoherent scatter spectrum, including the enhanced plasma lines assuming 
Maxwellian secondary electrons, resulting in the first radar measurement of its 
flux and a current carried by the thermal electrons. 



2 Experiment 



The measurements we present were collected by the 930 MHz EISCAT UHF in- 
coherent scatter radar, with its transmitter located at Ramfjordmocn in Norway 
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(69.6 °N, 19.2 °E, L=6.2). The signals scattered from the ionosphere were received 
at stations in Kiruna, Sweden and Sodankyla, Finland as well as at the transmit- 
ting site. General descriptions of the radar facility are given by Folkestad et al 



£1982) and Baron (1984 ). Local magnetic midnight at Ramfjordmoen is at about 



2130 UT. The aim, in the Swedish-Japanese EISCAT campaign in February 1999, 
was to measure the ionospheric parameters inside and outside the auroral arcs. 

For this a 3 channel ion line alternating code (Lchtincn and Haggstrom, 1987) 
experiment, optimised to probe the E-region and lower F-region with as high a 
speed as possible, was developed. The 16 bit strong condition alternating code with 
bitlengths of 22 fis was used, giving 3 km range resolution, and with a sample rate 
of 11 /is the range separations in consecutive spectra were 1.65 km. Fig. [l] shows 
the transmission/reception scheme of the first 20 ms of the radar cycle. The whole 
alternating code sequence takes about 0.3 s to complete. During this period the 
incoherent scatter autocorrelation functions (ACF) at the probed heights should 
not change significantly for the alternating codes to work. In order to keep this 
as short as possible, the short pulses, normally used for zerolag estimation, were 
dropped and instead a pseudo zero lag, obtained from decoding the power profiles 
of the different codes in the alternating code sequence, was used. Fig. || shows the 
range-lag ambiguity function for this lag centred around fj,s in the lag direction, 
but the main contribution to the signal comes from around 7 /is. The more normal 
lag centred at 22 fis is also shown for comparison. The range extents are rather 
similar but the power is, of course, considerably lower for the pseudo zero lag. 
Nevertheless, this is taken care of in the analysis and this lag is rather important 
in events with high temperatures giving broad ion line spectra or narrow ACFs. 
The transmitting frequencies were chosen to give maximum radiated power for a 
given high voltage setting. 
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Fig. 1. Transmission (black) and reception (dark gray) scheme for part of SP-SW-ALIS. A 
16-bit alternating code, 22 /is bits, is cycled over 3 frequencies. The interscan period is 9 ms and 
the total cycle takes 300 ms. The plasma line channels were set to sample the same range extent 
as the ion lines. 



The monostatic plasma line part of the experiment used two channels covering 
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Fig. 2. Range-lag ambiguity function for the first two lags, a) ion line and b) plasma line. The 
range, given in /is, can be converted to km with multiplication of 0.15. The differences between 
the ion and plasma line ambiguity functions are due to the use of different receiver filters. 



Channel setup for SP-SW-ALIS 990214 
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Fig. 3. The frequency setup for the experiment. The plasma line channels were fixed to 926 
and 934.5 MHz, giving "simultaneously" up and downshifted frequencies at 3 and 6.5 MHz. In 
addition there is an upshifted line at 4 and a downshifted at 5.5 MHz. 
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the same ranges as the ion line but with 3.3 km range separations between the 
spectra. The frequency setup for the experiment, illustrated in Fig. gives the 
possibility for 3 upshifted and 3 downshifted bands. For 3 MHz and 6.5 MHz 
frequency shifts, both the up- and downshifted plasma lines were measured. In 
addition there was a 4 MHz upshifted and a 5.5 MHz downshifted band. The 
width of these bands should have been set to match the bit length of the codes 
used, 50 kHz, but unfortunately this was not the case and 25 kHz wide filters 
were erroneously used. This gave naturally less signal throughput, and in Fig. || 
the corresponding range/lag ambiguity functions for the plasma line channels are 
included for the first two lags. The decoding still works, giving just slightly in- 
creased unwanted ambiguities, but above all the pseudo zero lag is moved out to a 
larger lag value. This fact, with one exception, almost ruled out the possibilities to 
measure plasma line spectra because, as will be shown here later, they are likely to 
be rather wide, due to large time and height gradients of the plasma parameters. 
This lead the analysis to use mainly the undecoded zerolag, which after integration 
over the different codes, almost resembles the shape of a 352 [is (16x22 [is) long 
pulse. 

The experiment contained a large number of antenna pointings in order to follow 
the auroral arcs, but as this day was cloudy over northern Scandinavia the trans- 
mitting antenna was kept fixed along the local geomagnetic field line. The remote 
sites, receiving only ion lines, were monitoring the same pulses as the transmitting 
site, and were used to measure the drifts in the F-region, to derive the electric field. 
Thus, these antennas intersected the transmitted beam at the F-region altitude 
giving the best signal, for this day mostly at 170 km. 



3 Measurements 

3.1 Ion line 

The experiment started at 1900 UT on 14 February 1999 and continued until 
2300 UT. Fig. |i| shows an overview of the parameters deduced from the ion line 
measurements, which were analysed using the on-line integration time, 5 seconds, 
in order to be compared to the plasma line data. The analysis was done using 



the GUISDAP package (Lchtincn and Huuskoncn, 1996), but a correction of 45% 



of the radar system constant used in the package had to be invoked to fit the 
plasma line measurements according to Eq. ([!]). This short integration time was 
possible due to the highly optimised mode used, with all the transmitter power 
concentrated to the E-region. In range, some integration was done, so that at 
lower heights 2 range gates were added together and with increasing height the 
number of gates added together increases to 15 in the F-region. 

There is a rather strong E-region from the start, but no real arcs, and we 
interpret this as diffuse aurora. The peak electron density shows some variation, 
but as time goes the E-region ionisation decreases until 2050 UT, where it is almost 
gone. The density peak during this time was at around 120 km altitude and the 
lower edge of the E-region at 110 km, but at times the ionisation reaches down to 
100 km. At 2050 UT and onwards until 2240 UT the ionosphere above Troms0 
became more active and several auroral arcs passed the beam. Around some of the 
arcs there are short-lived enhancements of electric field, seen as F-region ion and 
E-region electron temperature increases. In the last 10 minutes of the experiment 
the arc activity disappeared and again there was diffuse aurora. From the field 
aligned ion drifts it is evident that there is a rather strong wave activity in the 
diffuse aurora until 2050 UT, while it is not so clear in the continuation of the 
experiment. The last panel with the inferred electron density from the pseudo 
zero lag shows the same features as the fitted density panel but with highest 
possible resolution since no height integration is made. 
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Fig. 4. Summary of results from the ion line experiment. All panels shows the parameters 
in a altitude versus time fashion with the antenna directed along the geomagnetic field line. 
Panels from top: Electron density, electron temperature, ion temperature, line-of-sight ion drift 
velocity and electron density based on only the returned power. The data were analysed with 
5 s integration due to the active conditions. 
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3.2 Plasma line 



Since the analysis of the plasma lines was forced to handle the undecoded zerolags 
of the alternating codes, it was necessary to analyse their profile shape. Fig. || 
shows how this analysis was performed. From the fitted parameters of the ion line, 
electron density and temperature, a profile of the approximate Langmuir frequency 
can be calculated using Eq. (|l|). When probing at fixed frequency, there will be 
scattered signal only from heights where the probing and Langmuir frequencies 
match each other. The effect of the undecoded zerolag is similar to the one where 
a normal long pulse is used, but with a lower signal strength. So, the profile shape 
will be a square pulse centred on the corresponding altitude, since no gating is 
performed. Because the signal strengths are rather weak compared to the system 
temperature, there is a great deal of noise in the profile shape. To be able to 
extract the altitude and signal power, a fitting procedure need to be performed. 
For this a continous function and a good first guess is needed and the measured 
plasma line profile was convolved with the pulse shape to have a triangular shape 
and also showing a peak close to the matching height. 
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Fig. 5. Description of the analysis procedure, a) Langmuir frequency profile together with the 
probing frequency of 3 MHz shift, b) The echo profile using a short pulse, c) The echo profile 
using a long pulse, d) The echo profile of a long pulse convolved with the pulse shape. 



The plasma line part of the experiment is overviewed in Fig. |6J. The signal 
strengths shown should be compared with the UHF system temperature of about 
90 K. At first glance there is almost nothing in the upshifted part, but a more 
careful look shows weak signals between 1940 and 2030 UT and after 2250 UT, 
corresponding to the occurrence of diffuse aurora. Similar echoes can also be 
seen in the downshifted channel, and are due to plasma lines at 3 MHz offset 
from the transmitted frequencies, according to the ion line measurements. In the 
downshifted part after 2100 UT, frequent events of rather strong signals in phase 
with auroral arcs pass the beam. Most of these are from the 5.5 MHz shift, but 
some of them also are due to plasma lines at 6.5 MHz. Such events are less frequent 
and for most of them there are also signals in the upshifted part. 

Due to the fact that the same channel is used for several frequencies, there 
can be several altitudes that fulfill the matching condition between Langmuir and 
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Fig. 6. Overview of the plasma line measurements. The panels show the undecoded upshifted 
(top panel) and downshifted returns, given as signal strengths in Kelvin. 



probing frequency. This complicates the analysis somewhat, and there have to be 
several fits with different numbers of triangles superposed on each other. Some 
examples of this analysis are shown in Fig. ^, from only one plasma line to several 
both up- and downshifted lines. Using a lower limit of 2 K signal power, the total 
number of 5 second integration events with enhanced plasma lines for this evening 
was 256, and a total of 468 plasma line echoes were detected, divided into 220 on 
3 MHz, 19 on 4 MHz, 157 on 5.5 MHz and 72 on 6.5 MHz. 



4 Theory 

In order to relate the plasma line measurements to physical quantities it is neces- 
sary to investigate the spectrum of the incoherent scatter process. The Nyquist 
theorem approach, derived in a long series of papers by Dougherty and Farley 
1961 [L963D , |Farley et al. (196l| ), |Farley (196(1 and finally [Swartz and Farley 



(197G ), arrives at 
where 



N e r 2 sin 2 5 
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(2) 



N e e 2 



(3) 



and the index e stands for electrons and i for the different ion species. N and n 
are the densities, r e the classical electron radius, v the bulk velocity q the charge, 
e the electron charge and T the temperature. The complex normalised admittance 
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Fig. 7. Examples of measured convolved profiles (circles) and the triangle(s) fitted (lines). 
Included on the plots are the profiles of the Langmuir frequency (black) calculated from the 
parameters of ion lines fit. The value on the x-axis represent signal strength (K) and frequency 
(MHz) respectively. The numbers in the title shows the time of measurement and how many 
echoes was detected for the up- (red) and downshifted (green) frequencies. A lower limit of 2 K 
received signal power (2% SNR) was used, a) Simultaneous up and downshifted echoes at 3 MHz 
in diffuse aurora, b) An upshifted 3 MHz echo in diffuse aurora, c) A strong downshifted echo 
at 5.5 MHz in an auroral arc. d) Simultaneous up and downshifted echoes at 6.5 MHz and a 
downshifted 5.5 in an auroral arc. 
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function, y, contains most of the physics with the plasma dispersion function and 
have as main arguments th e collision frequency and magne tic field. The same 
result was also reached by Rosenbluth and Rostoker (1962 ), using the dressed 
particle approach. In Fig. || there is an example of the spectrum, showing clearly 
the strong ion line around zero offset frequency and the rather weak plasma lines 
at rather large offsets. The plasma lines become enhanced by a photo electron 
or auroral electron produced suprathermal electron distribution, and in order to 
simulate what this extra distribution does to the spectrum a modification to the 
formula has to be made. First, a rewriting of Eq. (||) following Swartz (197S| ) has 
to be performed in order to separate the electron and ion contributions to the 
spectrum: 



S(f) = - 

TT 



f+kvi/2-x 



Ti 



f+kv c /2ir 



where 



C D = 



k 2 eK 



niZfyi 



Ti 



(5) 



(6) 



/ is the frequency shift, Z — q/e, e is the dielectricity and K is the Boltzmann 
constant. Here, a normalisation of the spectrum has also been performed, so that 
the zero lag of the corresponding ACF is the raw electron density and the vector 
velocity is replaced by the line-of-sight velocity. Using a treatment in analogy 
to the ion contribution, it is now possible to rewrite the spectrum to support a 
number of Maxwellian electron distributions as 
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(7) 



In Fig. U the effect on the plasma lines of a suprathermal distribution with a 
reasonable density of 10 7 m -3 and width of 10 eV is shown, being lower than 
the ionisation energy for most ions. The plasma lines grow considerable and the 
integrated power over the bandwidth used in the experiment becomes comparable 
to the power of the ion line. It may also be noted that there is no effect at all seen in 



the ion line. Perkins and Salpeter (1965) have shown similar calculations, but when 



their method was based on large expansions to allow non-Maxwellian distributions 
one can here more directly superpose a few Maxwellian distributions to explain 
the measurements and even make fits of the spectra taken to get estimates of the 
suprathermal distributions, which is of great importance in auroral measurements. 
A consequence of Eq. (Q) is that it is also possible to derive the spectrum assuming 
currents carried by the suprathermal distribution, since it allows different drift 
velocities on the various distributions. Indeed, Fig. [l(], shows differential strengths 
on the two plasma lines, with upgoing electrons enhancing mainly the downshifted 
line and downgoing ones the upshifted line. 



5 Discussion 

Plasma line measurements in the active auroral ionosphere are not an easy task, 
due to the large variations in the ionospheric parameters. The Langmuir frequen- 
cies are largely dependent on the ambient electron density, making the line move 
considerably as the density changes, which it does on time scales of seconds. More- 
over, the density height gradient makes the lines very broad when measured over 
a specific height interval, and at times even broader than the receiver band. A 
chirped radar would solve only a part of the problem at the cost of transmitter 
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Fig. 8. The wide frequency incoherent scatter spectrum. The parameters used: N e = 4 ■ 
10 11 m -3 , T e = 700 K, Tj = 600 K, Vi = v e = ms" 1 and collision frequencies u e and Vi using 
the MSIS90e model (Hedin, 1991) for 140 km altitude. Logarithmic scale for the strength is used 
to be able to emphasise the different lines of the spectrum. 



power. These complications make it hard to draw any conclusions on the power in 
the lines, as one does not know for sure the scattering volume or the time duration 
of the scattering. Bearing this in mind and to at least minimise these effects, one 
can nevertheless look at the different distributions in height and power for the 
different lines to get an idea of their nature, using the on-line integration time of 
5 seconds. 

The altitude distributions for the different frequencies are shown in Fig. [n]. 
One must note that the signal levels shown are not corrected for range, as the 
scattering volumes are not known, so signals from a higher altitude are in fact 
stronger than the corresponding signal from a lower height. It is clearly seen 
that the strongest signal is the 5.5 MHz line, followed by 6.5, 4 and 3 MHz. 
The altitude distribution shows more or less the expected dependence on range, 
but there are some exceptions: In one point at 188 km in the 3 MHz band and 
for the 5.5 MHz band the strong values between 130 and 140 km seem to be 
stronger than the others, even taking into account the range effect. However, the 
number of points are too few to be used as evidence on altitude effects. Most 
of the echoes are coming from around 120-150 km altitude and Fig. |l2| shows 
a simulation of the expected strength of the plasma line for given background 
electron density. It shows a peak at around 5.5 MHz and this is also what the 
experiment shows. A more realistic suprathermal distribution will decrease the 
returned power for a number of frequencies and one should see the figure as an 
upper-limit estimate. Indeed, Nilsson ct al. (1996b| ), have made predictions of the 
expected strength of Langmuir waves for different heights and carefully derived 
distributions. These predictions are in rather good agreement with the present 
measurements showing a strong peak between 5 and 6.5 MHz. Although there is 
some uncertainty on the scattering volume, the observed strengths of the plasma 
lines can be used to get some estimates of the distribution of the suprathermal 
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Fig. 9. Close-up of the ion and the plasma lines. Top panel shows the electron distribution, 
middle panel the ion line and bottom panel the plasma lines, where the frequency scale of 
the downshifted line have been reversed. The figure shows the lines with two different electron 
distributions, thick line for the normal thermal distribution and thin line for the same distribution 
together with an suprathcrmal distribution. 
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Fig. 10. The plasma lines based on an upgoing (thick line) and downgoing 10 eV beam of 
suprathermal electrons at a current of 1 (iAm -2 and temperature of 5 eV. The thermal param- 
eters are the same as in Fig. H 



population. Assuming a width of 6 eV, one need, to get to the measured strengths 
for the 3 MHz suprathermal density of about 5 • 10 8 m 3 which is 0.5 % 

of the thermal population. For the 5.5 and 6.5 MHz cases, it is enough with only 
1 • 10 8 m -3 , but the uncertainty of volume is even larger here due to the active 
environment and the values should maybe be the same as for the lower frequency 
offsets. 

The most interesting thing with incoherent plasma lines is, of course, the pos- 
sibility to derive differential drifts between ions and electrons, and from these 
deduce ionospheric currents. For this, one needs to measure the up- and down- 
shifted lines simultaneously. Although the time of the measurement for both lines 
were not exactly the same in this experiment, the time shifts between them are so 
small (3-6 ms) compared to the total cycle time (300 ms) of the codes, that this 
effect is of minor importance. In Fig. |l3j, the strengths and altitudes of the two 
concurrently recorded up- and downshifted plasma lines are shown. In general, 
there are stronger up- than downshifted lines for the 3 MHz case, whereas no such 
trend can be seen in the 6.5 MHz band. However, there are exceptions to these 
overall trends and on occasions there are large differences in the signal strengths 
between the lines. Almost all of the 3 MHz plasma lines were recorded in dif- 
fuse aurora and this evident difference in signal power needs a closer examination. 
When there is a drift of the thermal electrons, the plasma lines shift, and when 
probing at a fixed frequency, the scatter may not come from the same altitude 
for the up- and downshifted lines respectively. The strength of the plasma line is 
also rather altitude dependent due to the damping by the collisions of electrons 
with ions and neutrals. But Fig. |l3| shows no general height difference between 
the up- and downshifted 3 MHz lines, so this difference in strength cannot be ex- 
plained by thermal electron bulk drifts. To simulate the effect of current carried 
by suprathermal electrons, Fig. [l4| illustrates the strength of the 3 MHz plasma 
lines for Maxwellian electron beams of different energies. With no net current the 



14 



3 MHz 



4 MHz 



5.5 MHz 



6.5 MHz 



50 




Power (K) 



Power (K) 



50 
Power (K) 



Power (K) 



Fig. 11. Scattered plots of echo altitude and strength for the different frequency shifts. The 
altitudes range from 100 to 210 km, E and lower F region. 
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Fig. 12. Simulation of the plasma line strength at 140 km altitude with an added distribution of 
secondaries for different electron densities from 0.5 to 10 ■ 10 11 m -3 . The envelope in the lowest 
panel shows the total power of the plasma lines versus the Langmuir frequency. 
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lines are of almost the same strength, and the difference is mainly dependent on 
where in the receiver band the lines are. However the experiment shows stronger 
upshifted lines, thus it is evident that the diffuse aurora this night contained fluxes 
of downgoing suprathermal electrons or, in other words, there was an upgoing cur- 
rent carried by suprathermal electrons. The average power ratio between up- and 
downshifted lines, 1.4, can be used to estimate the amount of current carried by 
the suprathermals. As before, using a density of 5 • 10 8 m~ 3 and temperature of 
6 eV, and shifting this population to a current density of 15 /xArn" 2 brings to the 
observed ratios. This current seems somewhat high, but not unrealistic. 

For the 6.5 MHz bands there may be a slight difference with respect to the 
altitude, and that is most likely due to thermal currents causing frequency shifts 
of the plasma lines. This effect is not very clear, but as lower heights have higher 
density, or Langmuir frequency and as the downshifted line is at a slightly lower 
height, this is most probably an effect of a downgoing current carried by thermal 
electrons. The correction due to heat-flow in the plasma dispersion function, 
discussed by Kofman et al. (1993) and confirmed later also by Nilsson et al. (1996a| ) 
and Guio et al. (1996), but not taken into account here, would also show the same 
effect in altitude difference between the lines. The band widths used here, 25 kHz, 
are much wider than the effect of heat-flow, which is less than or around 1 kHz in 
the F-region and much lower in the E-rcgion, so that c annot explain the 6.5 MHz 
height differences. A recent paper by Guio et al. (199§| ) with proper calculation of 
the dispersion equation investigates the heat-flow and finds that it is not necessary 
to invoke the effect at all, but as this paper don't go to the same extreme the heat- 
flow has to be considered. 
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Fig. 13. Scatter plots of power and altitude for the two simultaneously measured up and 
downshifted lines at 3 and 6.5 MHz. The triangles and circles for the power of the 3 MHz case, 
in the upper left, show lines detected before and after magnetic midnight respectively — no 
difference can be seen. 
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Fig. 14. Simulation of 3 MHz plasma line strength at 140 km with an electron beam with 
varying mean energy from to 20 eV and temperature of 10 eV, corresponding to a few /iAm~ 2 
field-aligned currents, depending on energy. The top panel shows the ratio between plasma and 
ion line strengths for upgoing beams and the bottom for downgoing beams. 



Most of the above discussions on currents are only on directions, but to get 
any quantitative numbers it is necessary to look at the individual spectra. Of 
the 256 events of enhanced plasma lines, there is only one that is good enough 
to investigate. All the other are too broad either due to the Langmuir frequency 
gradient smearing or the fact that the Langmuir frequency at the height in question 
is varying during the 5 second time slots. A further reason is that the stationarity 
condition for the alternating code technique is not fulfilled, due to changes of 
Langmuir frequency, and the spectra change too much within the 300 ms cycle. 
Anyway, there is one, and the up- and downshifted bands, together with the ion 
line band, are shown in Fig. |l5. The shifts in this case are around 3 MHz as 
deduced from the ion line analysis. These spectra were then fitted to the theoretical 
spectrum in Eq. (R) using 7 ionospheric parameters, the thermal parameters for 
both ions and electrons and the parameters of a non-shifted suprathermal electron 
distribution. At the first glance the fit on the ion line seems rather poor, but the 
fit was done as usual in the time domain with proper weighting on the different 
lags of the ACF, and in the FFT process to produce Fig. [jj] these statistical 
properties are lost. The fit looks actually better in the time domain, but the 
frequency domain was chosen in the figure to be more informative. The current 
density j — N e e(vi — v e ), and with the fitted parameters, the field aligned current 
carried by the thermal electrons, amount to 12 fiAm" 2 downward. This is opposite 
the general current seen in the plasma line strength for the 3 MHz band, but 
the magnitudes are comparable. For this case no current could be seen on the 
suprathermals, and it is evident that the currents are rather structured. Of course, 
there are many more cases of no plasma lines than enhanced plasma lines in the 
diffuse aurora, and some cases of only one plasma line either up- or downshifted. 
It is not surprising that this single example of thermal currents do not follow 
the general trend. Again, heat-flow is not taken into account, but that should 
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Fig. 15. Measured spectra of the both plasma lines and the ion line and the best 7-parameter fit 
of the theoretical spectrum. The parameters fitted were N e = 9.7 ■ 10 10 m" 3 , T e = 631 K, T; = 
697 K, Vi = —100 ms^ 1 , v e = 657 ms" 1 and a suprathermal distribution with n e = 0.8- 10 9 m -3 
and T = 11 eV. 



18 



increase the current s omewhat. The deduced suprathcrmal distribution is in good 
agreement with those Kirkwood et al. (1995| ) derived from the precipitating flux 
of primaries; of course in the present case the distribution is Maxwellian and 
without the fine structure due to the atmospheric constituents, but the numbers 
are comparable. 



6 Conclusions 

We have made measurements of plasma lines in the active auroral E-region. During 
256 periods of 5 second integration we found a total of 468 plasma line echoes, 
divided into 220 on 3 MHz, 19 on 4 MHz, 157 on 5.5 MHz and 72 on 6.5 MHz. It 
may seem strange to try to measure plasma lines at such a low frequency as 3 MHz 
giving low signal levels, but in fact most of the echoes and the most interesting 
results came from this frequency offset. 

The strongest echoes were found at the 5.5 MHz line, and somewhat weaker 
ones at 6.5 MHz, inside auroral arcs. One must, however, note that the strength 
measured inside the arcs is mostly a low-limit estimate due to the active environ- 
ment. The integration period used, 5 s, is rather long in auroral arc conditions, 
and changes typically occur on shorter time scales. Therefore, effects of gradients 
in the Langmuir frequency profile, and hence scattering volumes, have not been 
taken into account. 

The Holy Grail in incoherent scatter plasma lines is the possibility to measure 
currents, and in this case, the field-aligned currents in aurora. The simulations car- 
ried out here, the extended full incoherent scatter spectrum with multi-Maxwellian 
distributions of electrons, show that the strength of the lines is determined by the 
suprathcrmal part of the electron distribution, and the frequency mainly by the 
thermal part. For simultaneous up- and downshifted plasma line differences in 
intensity we can deduce currents carried by suprathcrmals, and for differences in 
frequency, currents carried by thermals. 

The simultaneous up- and downshifted frequencies of the 3 MHz line in the 
diffuse aurora show, on average, an upward field-aligned suprathermal current 
during the two main periods when they were detected, 1940-2030 UT and 2250- 
2300 UT. In the arcs in general, there is an indication of downward thermal current 
as seen from the altitudes of the 6.5 MHz echoes. Of course, no rule is without 
exceptions, and there are cases where one line is much stronger than the other or 
the other line is not at all enhanced, indicating strong currents. 

In the full 7-parameter fit of the incoherent scatter spectrum with the ion 
line and the both enhanced plasma lines, we obtained a thermal current with 
a suprathermal distribution of electrons consistent with distributions derived from 
precipitating fluxes. 
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